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Improved Analysis of J/V' Decays into a Vector Meson and Two Pseudoscalars 
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Recently, the BES collaboration has published an extensive partial wave analysis of experimental 
data on J/tp <j}n^n~ , J/ip —> wtt+tt", J/ip cj>K^ K~ and J/t/j uK^K~ . These new results 
are analyzed here, with full account of detection efficiencies, in the framework of a chiral unitary 
description with coupled-channel final state interactions between tttt and KK pairs. The emission of 
a dimeson pair is described in terms of the strange and nonstrange scalar form factors of the pion and 
the kaon, which include the final state interaction and are constrained by unitarity and by matching 
to the next-to-leading-order chiral expressions. This procedure allows for a calculation of the S-wave 
component of the dimeson spectrum including the /o(980) resonance, and for an estimation of the 
low-energy constants of Chiral Perturbation Theory, in particular the large A'^c suppressed constants 
L\ and Lg- The decays in question are also sensitive to physics associated with OZI violation in 
the 0^"'' channel. It is found that the S-wave contributions to ^tt^tt", <j)K'^ K~ and uj-k^-k~ given 
by the BES partial-wave analysis may be very well fitted up to a dimeson center-of-mass energy 
of ~ 1.2 GeV, for a large and positive value of L\ and a value of ig compatible with zero. An 
accurate determination of the amount of OZI violation in the J/'4> ^ <j)n^Ti~ decay is achieved, and 
the S-wave contribution to luK^ K~ near threshold is predicted. 

PACS numbers: 13.20.Gd, 12.39.Fe 
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I. INTRODUCTION 



The decays of the J /tp into a vector meson such as 
(j) ov Lo, via emission of a pair of light pseudoscalar 
mesons, may yield insight into the dynamics of the 
pseudo-Goldstone bosons of QCD 0, 0, 0] , and in par- 
ticular into the final state interaction (FSI) between tttt 
and KK pairs, which is an essential component in a 
realistic description of the scalar form factors (FFs) of 
pions and kaons. Additionally, such decays can yield 
insight into violation of the Okubo-Zweig-Iizuka (OZI) 
rule H H in the scalar (0++) sector of QCD, smce 
the leading order contributions to such decays are OZI 
suppressed. Furthermore, as shown in Fig. ^ a doubly 
OZI-violating component may contribute to these J /jb 
decays, which was demonstrated already in Re fs. B, Q 
although the data from the DM2 0, MARK-III |ld| and 
BES 1 1] collaborations available at that time had rather 
low statistics. However, since then the BES collabora- 
tion has published far superior data on J /ij) — > ^tt+tt" 
and J/ip (j}K'^K~ lial. as well as for J/i/; — > wtt+tt^ 
and J/V' ^ wK+fsT- [ll[l3. Additionally, a compre- 
hensive partial- wave analysis (PWA) has been performed 
for those data, which is particularly significant since an 
explicit determination of the S'-wave tttt and KK event 
distributions is thus available. Thus, a much more pre- 
cise analysis of the issues first touched upon in Ref. 
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clearly called for. 

A key ingredient in such an analysis is a realistic treat- 
ment of the final state interaction (FSI). It has been 
demonstrated in Ref. that the FSI in the tttt — KK 
system can be well described by a coupled-channel Bethc- 
Salpeter approach using the lowest order CHPT ampli- 
tudes for meson-meson scattering In such an 
approach, the lowest resonances in the 0"*""*" sector are of 
dynamical origin, i.e. they arise due to the strong rescat- 
tering effects in the tttt or KK system. Such dynamically 
generated states include the a and /o(980) mesons, which 
are prominent in the BES data on dimeson emission from 
the J/V' [Illlllll. It is useful, in view of the controver- 
sial nature of the /o(980) 0, H S IH US HI ilh to recall 
Ref. , which generalizes the work of Ref. |l5l| . There, 
explicit S'-wave resonance exchanges were included to- 
gether with the lowest order CHPT contributions in a 
study of the partial wave amplitudes for the whole scalar 
sector with 1 — 0,1/2 and 1. It was noted that the results 
of Ref. 15] could be recovered when the explicit tree- level 
resonance contributions were dropped. The conclusion of 
Ref. was that the lowest nonet of scalar resonances, 
which includes the a, k and ao(980) states, is of dynami- 
cal origin, while a preexisting octet of scalar resonances is 
present at ~ 1.4 GeV. It was also noted that the physical 
/o(980) state obtains a strong contribution of dynami- 
cal origin, and may also receive one from a preexisting 
singlet state. 

This analysis uses the formalism introduced in Ref. j3| , 
where the expressions for the scalar FFs of the pions and 
kaons were obtained using the results of Ref. _15j. This 
allows for a description of the scalar FFs which takes 
into account the FSI between pions and kaons up to 
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II. THEORETICAL FRAMEWORK 

The theoretical tools required for the calculation of the 
scalar FFs of the pseudo-Goldstone bosons using CHPT 
and unitarity constraints have, as discussed in the Intro- 
duction, already been extensively treated in the existing 
literature, and therefore only the parts directly relevant 
to the present analysis are repeated here. For conve- 
nience, the NLO expressions for the scalar FFs of the 
pseudo-Goldstone bosons are explicitly given here. Also, 
the expressions for the scalar FFs of the pion in Ref. 
require a minor correction^, such that an updated version 
is called for. A derivation of the scalar FFs using unitar- 
ity and the methods of Ref. [l^ is given in Ref. , and 
introductions to CHPT can e.g. be found in Ref. 



A. Amplitude for tttt and KK emission 



FIG. 1: Quark line diagrams for the decay of a J/ip into 
a light vector meson (0) and a pair of light pseudoscalar 
mesons {tt'^-k~). The annihilation of the cc pair proceeds 
via at least three gluons. The shaded squares denote the 
strange (F^) and nonstrange (F^) scalar FFs of the pion. Note 
that diagram (B) is further OZI-suppressed with respect to 
diagram (A). 



~ 1.2 GeV. At higher energies, a number of preexist- 
ing scalar resonances such as the /o(1500) have to be 
accounted for, as well as the effects of multiparticle in- 
termediate states, most importantly the Att state. These 
scalar FFs may then be constrained by matching to the 
next-to-leading-order (NLO) chiral expressions. This al- 
lows for a fit of the large Nc suppressed Low-Energy Con- 
stants (LECs) L\ and Lg of CHPT to the dimeson spec- 
tra of the J/ip ^ (pTTTT and JM (pKK decays, using 
the Lagrangian model of Ref. Q ■ The amount of direct 
OZI violation present in these decays may also be accu- 
rately estimated. Finally, it should be noted that the 
present treatment of the FSI has been proven successful 
in describing, in the same spirit, a wide variety of pro- 
cesses, such as the photon fusion reactions 77 — > tttt and 
77 — !■ KK p?! , the decays — > 77r7r and 4> ^ ^KK , 
and more recently the hadronic decays D, Ds —>■ tttttt and 
D,Ds Ktttt 26j. 

This paper is organized in the following manner: In 
Section ^ the description of the J/tp (piriT decays is 
briefly reviewed, along with the FSI in the tttt — KK sys- 
tem, as applied to the scalar FFs of the pseudo-Goldstone 
bosons. Some minor corrections to the NLO CHPT ex- 
pressions in Ref. are also pointed out. Section [ffl] de- 
scribes the analysis of the experimental BES data, along 
with a discussion of the fitted parameter values, with 
emphasis on the LECs of CHPT and the evidence for 
OZI violation. In Section Hvl the results are summarized 
along with a concluding discussion. 



This work makes use of the SU(3) and Lorentz invari- 
ant Lagrangian of Ref. Y] to describe the decay of a J/ip 
into a pair of pseudoscalar mesons and a light vector me- 
son. This Lagrangian can be written as 



C = 



i.Vf El) , 



(1) 



where the Vs and Vi denote the lowest octet and singlet of 
vector meson resonances. Similarly, the Ss and Ei refer 
to the corresponding sets of scalar sources, as defined in 
Ref. 0. In the above equation, g denotes a coupling 
constant, the precise value of which is not required in 
the present analysis, while the real parameter ly will be 
shown to play the role of an OZI violation parameter in 



the 



channel. Furthermore, the angled brackets 



in Eq. ^ denote the trace with respect to the SU(3) 
indices of the matrices Vs and Sg. Evaluation of that 
trace yields 



(2) 



where only the terms relevant for the present considera- 
tions have been written out. Here Vg denotes the 1 = 
state of the octet of vector meson resonances, while Ss 
again refers to the corresponding operator in the matrix 
of scalar sources. The (f> and uj fields, along with the as- 
sociated scalar sources Scf, and S^j are then introduced 
according to 



f2 



V3' 



2 c. 



(3) 
(4) 



^ The scalar FFs of the pion in Ref. 0| require a minor correc- 
tion in the values at s = 0. The authors thank J. Bijnens and 
J. A. Oiler for their assistance in pinpointing this. The effect on 
the numerics of Ref. is negligible. 
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which corresponds to the ansatz of ideal mixing between 
Vg and Vi- The departure from this situation is reviewed, 
using different models, in Rcf. 28]. The amount of de- 
viation from ideal mixing in the cf) — uj system has been 
estimated (21 to influence, in an analysis of the present 
kind, the determination of the magnitude of the OZI vi- 
olation at the ~ 5% level. This should be compared with 
the expected ~ 40% departure from unity of the pa- 
rameter V in Eq. (01. In view of this, the relations given 
in Eqs. Q and Q will be considered adequate for the 
present analysis. 

The scalar sources and Suj may, in terms of a 
quark model description, be expressed as = ss and 
= fin = [uu + dd)/^/2. By means of these relations 
and Eqs. and the Lagrangian of Eq. ((2Jl may be 
written in the form 

+ ^,,u:^^C^{v)[ss + \^{v)nn], (5) 

where the coupling constant g is taken, as further elabo- 
rated in Sect. Ill B1 to be absorbed into and C^j- These 
and the \i in Eq. (O are given in terms of the parameter 
v according to 



_ V2{v - 1) 



V 

1 + 2v 
V2(i^-1)' 



V2(i. - 1) 



(6) 



(7) 



which shows that the parameters for the lo decay operator 
can be expressed in terms of those which control the </) 
decay. The explicit relations are given by 



A.., = 



(8) 



From now on, the dependence of the Ci and \i on v will 
be suppressed. The quantities to be determined from fits 
to the experimental dimeson spectra of Refs. [l^lTsLIT^ 
are taken to be and A^. It is worth noting that the 
limit A^ = corresponds to the value v ~ \, and in 
that case the dimeson spectra for </) and lo decays are 
driven entirely by the strange scalar source ss and the 
nonstrange scalar source nn, respectively. 

From the Lagrangian in Eq. |(SJ), the matrix elements 
and (pR'^K" decay of the J/tp are given by 



MY = A/;C^(0|(ss-f A^nn)K^)i*^o, (9) 
C^{0\{ss + X^nn)\KK)U>, (10) 



in terms of the tttt and KK states with / = 0, which are 
related to the physical tt+tt" and K~^K~ states by the 
Clebsch-Gordan (CG) coefficients in front of the above 
expressions. It should be noted that in Ref. 0, the 
CG-coefficient for tt+tt" decay was incorrectly written 



M 



M 



T 



FIG. 2: Overview of the decay J/tp <f)Tx^ii~ . The open 
square represents a vertex of the Lagrangian in Eq. Q, 
whereas the solid square denotes the scalar FF of the pseu- 
doscalar meson pair MM to NLO in CHPT. The shaded 
square represents the T-matrix for meson-meson scattering, 
as calculated from the LO CHPT amplitudes via the Bethe- 
Salpeter equation. The interpolating field S corresponds to 
the strange and nonstrange scalar sources ss and nn, as de- 
scribed in the text. 



as ■\/4/3. The full transition amplitudes also contain 
the polarization vectors of the J /ip and (p mesons, which 
have not been included in the above definitions. They 
introduce an additional, weakly energy dependent fac- 
tor which is given explicitly in Sect. ITT^ The matrix 
elements for cjtt+tt" and toK^K^ , may be obtained by 
replacement of the labels in Eqs. © and (|10|l according 
to (p LO. The matrix elements of the scalar sources are 
given in terms of the 7 = scalar EEs which are discussed 

in sect.rrrci 



B. Decay rates and dimeson event distributions 

The differential decay rate of a J/ip into a vector meson 
and a pair of pseudoscalar mesons is, for the case of ^tttt 
decay, given by 



dT 



(27r)3 



i^poi \M 



TTTT I 2 



(11) 



where PFn-Tr — \fs^. The decay rates for the other com- 
binations of (p and LO mesons and ixix^KK final states can 
be obtained by appropriate replacement of the indices in 
Eq. . The moduli of the (p and tt momenta in Eq. Hll|) 
are given by 



\V4>\ 



El 



2VF. 



(12) 
(13) 
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FIG. 3: Detection efficiency rj from BES Monte Cario simulations for J/^p <jf>7r"'"7r~ , J/ip <f)K^ K~ , J/ip u!-r^-r~ and 
J/ij) — > LoK^ K~ , along witfi fitted polynomial approximations for each case. The estimated detection efficiencies take into 
account the characteristics of the BES detector as well as the different cuts applied in the data analysis. The BES data 
published in Refs. IT^ are uncorrected for efficiency, and for that purpose the efficiency functions rj have been included 
in Eq. ifT^ . 



in the rest frame of the tttt system. The factor Fpoi in 
Eq. l(TT)) depends on the dipion energy Wt^-k and is gener- 
ated by properly averaging and summing over the polar- 
izations of the J /ip and </> mesons, respectively. It may 
be expressed as 

^poi ^ \Y.^,{p)e^{p)e:{p)e^*{p') 



Mi - Wi 



(14) 



Again, it should be noted that the corresponding expres- 
sions for Fpoi for the other decay channels considered can 
be obtained straightforwardly from Eq. 114|) by the sub- 
stitutions — > oj and TTTT — > KK. 

The results for Jjtp decay into a vector meson and 
a dimeson pair published by the BES collaboration are 
given in terms of event distributions as a function of 
the dimeson center-of-mass energy W . The relation be- 
tween the differential event distribution and the decay 



rate given by Eq. is defined to be 



dW^ 



dW^ 



(15) 



where the function ri{W) represents the detection effi- 
ciency, shown in Fig. |21 which also takes into account the 
effects of the various cuts imposed on the data in order to 
reduce the background to an acceptable level. It should 
be noted that the detection efficiencies ij cannot be ne- 
glected, since it is evident from Fig. O that a sizeable 
difference exists between the efficiencies for all four de- 
cays considered. Furthermore, the detection efficiencies 
exhibit significant nonlinear behavior. 

The overall constant of proportionality in Eq. (|15|l is 
not relevant to the present analysis, since it may be ab- 
sorbed in the definitions of and Ci^, along with the 
coupling constant g of Eq. Q and a factor V2Bo from the 
scalar FFs in Eq. (|19|l . Thus, the quantities with ^/2gBo 
and the proportionality factor absorbed are denoted by 

and While and gBo are dimensionless, 
has dimension [E^2^, 
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C. Scalar Form Factors from CHPT 

The nonstrange and strange scalar FFs of the pseudo- 
Goldstone bosons of CHPT are defined in terms of the 
S'-wave states with / = 0, 



/i=o 



lw)i=o 



i=|.+,-) + _L|,o,o^, (16) 
^\K^K-) + l=\K^KO), (17) 



71 '""^ 



(18) 



where |7r"'"7r ) denotes the symmetrized combination 
of Itt+tt") andj7r~7r+). Following the conventions of 
Refs. [s^, 1^, an extra factor of l/-\/2 has been in- 
cluded for the / = states composed of members of 
the same isospin multiplet. This takes conveniently into 



account the fact that the pions behave as identical par- 
ticles in the isospin basis. In terms of the above states, 
the scalar FFs for the 7r,K and rj mesons are defined as 



V2Bor^(.s) = 



(0\fm\i]?])^^Q, 



(19) 



where the notation {tt — 1, K — 2, rj = 3) has been 
introduced for simplicity. The expressions for the strange 
scalar FFs may be obtained by the substitutions F" — > F^ 
and nn —^ ss. As discussed above, the expressions given 
in Ref. are updated here with minor corrections to 
F" and Ff. With these definitions, the scalar FFs may 
be expressed in terms of the meson loop function Jfj 
and the tadpole factor /i^, given in Eqs. (|31|l and 
respectively. The expressions so obtained up to NLO in 
CHPT are, in agreement with Refs. [Ullllll, 



F?(s) 



for the pion, and 

ri(s) 



2 



l + M.-^ + ^(2XS-LD + 8(2Lg-LD^^^^^±M 



1 

71 



l+^(2.~m^-60 + ^(.-4mi) 



(6r 



' 9s - 8to|- \ jr I ^ 3s 3s 



8K 



riri\ 
2 /l™2 



8s 



4s 



16LE 



for the kaon. Finally, for the 77 one finds 



327^0 o 2 



(20) 
(21) 



(22) 



(23) 



ri(s) 



1 



1 



2V3 
128i 



24L^ 
P 



1 



10 



s + -m^ — 



3 " 
32L 



16 



S + -TO^ — 



4) 



9s — 8m|.\ 3 



6LI 
P 



16 



+ ^jf- - ml) + ^ {2ml - ml) - + 2/^ A' + 



4L£ / 4 2 16 2 \ 8Lg 2 2\ 
JTh + -ml - + -j^ {Sml + ml) 

16to|- - 7ml \ 

Igp J •^vny^' 



9g - 

) •JkkKs) 



(24) 



(25) 
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D. Matching of FSI to NLO CHPT 



where 



The constraints imposed by unitarity on the pion and 
kaon scalar FFs, the inclusion of the FSI via resum- 
mation in terms of the Bethe-Salpeter (BS) equation, 
the channel coupling between the tttt and KK systems, 
and the matching of the scalar FFs to the NLO CHPT 
expressions have all been elaborated in great detail in 
Refs. 0,0, and will thus be only briefly touched upon 
here. Within that framework, consideration of the uni- 
tarity constraints yields a scalar FF in terms of the alge- 
braic coupled-channel equation 



T{s) - [I + K{s)g{s)]-'R{s) 

= [I-K{s)g{s)]R{s) + 0{p% 



where in the second line, the equation has been expanded 
up to NLO, the NNLO contribution defined as being of 
0{p^) in the chiral expansion. This expansion is instruc- 
tive since it allows for the Jfj integrals from the NLO 
scalar FF expressions to be absorbed into the above equa- 
tion. Here K{s) denotes the kernel of S'-wave projected 
7 = meson-meson scattering amplitudes from the lead- 
ing order chiral Lagrangian. Using the notation defined 
in Sect. Ill CI they are given by 



2s -ml r. Vis 



K22 = 



2/2 
3s 



4/2' 



(27) 



where the constant / is taken to equal the pion decay 
constant /^, with the convention /^ = 0.0924 GeV. The 
components given above are sufficient for the two-channel 
formalism of Ref. used in this paper, where only the 
FSI in the tttt and KK channels is considered. The chiral 
logarithms associated with the rjr] channel can thus not 
be reproduced by Eq. (|26|l and are therefore removed 
from the chiral expressions given in Sect. Ill CI while the 
contribution of that channel to the values of the form 
factors at s = is retained. For completeness, it should 
be noted that if the rjr] channel is also included, then the 
matrix K(s) should be augmented by the elements 



Ki3 = K31 = 



K 



23 



K 



32 



2V3/2 
9s - 8m^ 
12/2 ■ 



i^.33 — 



7mi 



18/2 



(28) 



The elements of the diagonal matrix gi{s) in Eq. (|26|l are 
given by the cutoff-regularized loop integral 



9iis) 




CTiis) 



1 



imj 



(30) 



and ^max denotes a three-momentum cutoff, which has 
to be treated as an a priori unknown model parameter, 
which is however expected to be of the order of ~ 1 GeV. 
Since the above expressions are calculated in a cutoff- 
regularization scheme, it is useful to note that within the 
modified MS subtraction scheme commonly employed in 
CHPT, the meson loop function gi{s) is given by 



(26) J;;(s) 



1 



167r2 



1 - log 



a,{s) log 



<7i{s) + 1 

(Ti(s) - 1 



(31) 



for which it has been shown in App. 2 of Ref. |23 that 
an optimal matching between the two renormalization 
schemes requires that 



2gn 



(32) 



in which case the differences between the two forms are of 
0{m^/q^^^). Furthermore, the expressions for the log- 
arithms generated by the chiral tadpoles in the NLO 
scalar FFs are given by 



327r2/' 



(33) 



As demonstrated in Ref. |j7|, the quantity R{s) in 
Eq. (|26(l is a vector of functions free of any cuts or sin- 
gularities, since the right-hand or unitarity cut has been 
removed by construction. The information provided by 
CHPT can then be built into the formalism by fixing 
R{s) to the NLO CHPT expressions for the scalar FFs. 
Consideration of Eq. (|26|) yields the defining relations 



rr(s) = K(s) 



2 

J2K^,{s)g,{s)R;^{s)+0{p''), (34) 



where it is understood that only contributions up to 
0{p'^) are to be retained in the product KgR. The analo- 
gous expressions for the vectors Rf (s) associated with the 
strange scalar FFs rf(s) can be obtained from the above 
relations by the substitutions Ff Tf and i?f Rf. 
The above procedure is equivalent to the intuitive result 
obtained by dropping, in the expressions for the Fi, all 
occurrences of the loop integrals JJ^ and J^k^ and keep- 
ing only the parts of the J^^ which do not depend on s. 
Nevertheless, the explicit evaluation of Eq. H34I) provides 
a useful check on the consistency of the normalization 
used for the NLO scalar FFs and the LO meson-meson 
interaction kernel K. It should also be noted that the 
expressions for i?" and Rf correspond to the corrected 
scalar FFs, as explained in the beginning of Sect.|n] The 
expressions for the R^ so obtained are 
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3ml 



Rf{s) 



2887r2/2 
\/3 f Wml 



1 + log ^ 



P 



{2LI 



LI) 



8s 



727r2/2 




(35) 
(36) 



for the pion, and for the kaon one finds 




K 



K 



{iniK + ml) + 



32^70 2 2 



(37) 



(38) 



The expressions for the i?f and Rf are valid when only 
the TTTT and ifA' channels are considered in the FSI. On 
the other hand, if the full three-channel interaction ker- 
nel of Eq. (|27|) is used, then the above equations should 
be modified such that the terms in square brackets are 
dropped. With respect to the omission of the r/77 chan- 
nel, it was noted in Ref. that reproduction of the 
data on the inelastic tttt KK cross section requires 
the addition of a preexisting contribution to the /o(980) 
if the 7777 channel is included. On the other hand, no 
such contribution was found to be necessary if the 7777 
channel is dropped. Furthermore, the effect of this chan- 
nel is known |23| to be very small for energies less than 
~ 0.8 GeV. 

It should be stressed here, with respect to the above 
mentioned issues, that the main concern in the present 
analysis is the use of a meson-meson interaction kernel 
which is known to give a realistic description of the tttt 
phase shift close to the KK threshold. Since none of 
the adjustable model parameters have any influence on 
the behavior of the tttt phase shift, any model which has 
the proper chiral behavior for low energies and faithfully 
reproduces the /o(980) should therefore give similar re- 
sults. In view of this, the inclusion or omission of the 
7777 channel, or the question of a preexisting contribu- 
tion to the /o(980) resonance, are issues of secondary 
importance. Nevertheless, the uncertainties introduced 
by these issues into the determination of the L\ are dis- 
cussed in Sect. IIIII Finally, it should be noted that in 
order to minimize the dependence on m,,, the Gell-Mann 
- Okubo (GO) relation has been applied throughout in 
the polynomial terms of the and the Ri. 



III. FITS TO BBS DATA 

The event distributions given by Eq. (|15|l can be simul- 
taneously fitted to the dimeson spectra in the ^tt+tt", 
WTT+TT" and (pK^K^ channels. The parameters to be 
determined via a fit are the LECs L\,LI,Lq and Lg 
which influence the scalar FFs, as well as the model pa- 
rameters ^rf, and A^. Due to the accuracy of the BES 
data [12. 113L [Ijl , all of the model parameters can be well 
constrained, which is especially true for L4 and , while 
the sensitivity of the fit to L5 and Lg turns out to be 
somewhat less. Once all the model parameters are de- 
termined by a fit to the three decay channels mentioned 
above, the event distribution in the remaining loK^ K~ 
channel in essentially fixed. This channel is thus not in- 
cluded in the fit and is treated instead as a prediction. 
To a large extent, this also turns out to be true for the 
shape of the fitted unT'^-K~ distribution. 

In spite of the above mentioned positive issues, the fit- 
ting of the predicted event distribution dN/dW to the 
S-w&ve contribution from the PWA of the BES collabo- 
ration is complicated by several issues: Firstly, the de- 
tection efficiencies, determined by BES via Monte Carlo 
simulation and shown in Fig. |2| are different for each de- 
cay channel, and furthermore vary appreciably over the 
range of dimeson energies considered. Secondly, the S*- 
wave contribution in the BES PWA cannot be considered 
as strict experimental information, since it is inevitably 
biased to some extent by the parameterizations chosen 
there for the a pole. Thirdly, an unambiguous fit requires 
a highly precise description of the /o(980) resonance gen- 
erated by the FSI, which to a large extent cannot be ad- 
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justed in the present model. All of these issues, as well 
as the fitted parameter values and the associated error 
analysis are elaborated in detail below. 

A. Definitions 

In order for the fit results to be well reproducible, the 
various constant parameters which enter the expressions 
for the decay amplitudes should be accurately defined. 
These parameters include the masses of the light pseu- 
doscalar mesons, for which the current experimental val- 
ues ^3 have been used. These are = 0.13957 GeV, 
ruK = 0.49368 GeV and m,, = 0.5478 GeV, while the 
value /tt = 0.0924 GeV has been adopted for the pion 
decay constant. The physical masses of the charged pi- 
ous and kaons have been used in order for the tt+tt^ and 
K~^K~ thresholds to coincide with the physical ones. 
Also, the physical i] meson mass has been used rather 
than the one given by the GO relation. The t] meson 
mass appears in relatively few places in the expressions, 
and checks on the fits have indicated that replacement 
of the physical 77 mass with that given by the GO rela- 
tion has a minimal effect. Further parameters are the 
masses of the vector mesons p, w, and J/ip- The p mass 
is used in the evolution of the L^, and has been taken 
as nip = 0.776 GeV. The other vector meson masses 
appear in various phase-space factors, and have been 
given the values = 0.783 GeV, = 1.020 GeV 
and mj/^ = 3.097 GeV ia4i. 

It is not a priori obvious how the individual devia- 
tions, required for the ^t, are to be treated for the 
BES data. The individual deviations for each bin di, are 
given by the BES collaboration as the square root of the 
number of events These numbers represent the sta- 

tistical errors of the raw data, uncorrected for detection 
efficiency and background. Furthermore, what is fitted 
in the present analysis is not the total signal detected by 
the BES experiment, but rather the S'-wave contribution 
from the accompanying PWA. In view of these consid- 
erations, the deviations for each bin used in the fitting 
procedure have been taken to be of the form 

A, = V^^w,, (39) 

where the Wi represent weighting factors which have been 
chosen in a physically motivated way. In principle, indi- 
vidual Wi could be introduced in all the decay channels 
studied and would then represent a "quality factor" for 
each bin. In practice, to avoid excessive fine-tuning of the 
fit, a constant value for Wi has been applied for each de- 
cay channel, according to the following principles: Since 
the S'-wave contribution in the wtt+tt" spectrum is likely 
to have the largest uncertainty, a value of Wi = 3 has 
been adopted for that decay channel, whereas the val- 
ues of Wi for the 07r"^7r~ and (pK'^ K~ distributions have 
been defined to be between and 1. The choice of such 
a large value of Wi for wtt'^'tt" means that the shape of 
the 7r"'"7r~ spectrum is constrained by the <^7r"'"7r~ and 



(f)K'^K~ distributions, and may therefore be considered 
as a prediction of the present analysis. On the other 
hand, the relative magnitude of the ^tt+tt^ and lo'k^'k^ 
spectra provides an important piece of experimental in- 
put, without which all parameters in the fit cannot be 
constrained. Checks have been performed to verify that 
the choice of = 3 for u}'k^t:~ does not introduce un- 
physical skewing of the fit, and it was found that a further 
increase in wi for wtt+tt^ has a negligible effect on the fit 
results. The values for Wi used for the other two chan- 
nels are Wi = 1 for (/)7r+7r^ and Wi — 0.1 for (f)K~^K~ 
in "Fit I" and vice versa for "Fit 11" . Again, the fit re- 
sults have been checked for stability with respect to these 
choices. 

The predicted number of events Ni in each bin is fitted 
to the total number of S'-wave events determined by the 
BES PWA. To this end, the definition 

NiW) ^ J^) AW, (40) 

has been implemented, where the differential event distri- 
bution dN/dW is given by Eq. ifTs)) , and AW denotes the 
binsize chosen by BES for the PWA, given as 25 MeV for 
uj final states and 30 MeV for final states [ll[ll[l|. It 
should be reemphasized that all effects due to detection 
efficiencies are included in dN/dW as defined in Eq. lfTK|l . 
The number of events Ni per bin are then calculated 
by averaging Eq. l|in|) over the energy range [Wl,W2] 
spanned by the bin, giving 

1 r^i _ 
N. = ^l^NiW)dW, (41) 

where the explicit dependence on the binsize AW can- 
cels. The advantage of the function N{W) is that it can 
be conveniently compared with the experimental data, 
regardless of the binsize chosen for a particular decay 
channel. 

Finally, the determined in the fit are always quoted 
at the scale fj, = nip. The fit function itself depends on 
the LI at the scale given by Eq. H32|l in terms of the cutoff 
parameter Qmax of the FSI, which gives the L[ at a scale 
of /i ~ 1.2(3'max- For calculational purposes, the are 
evolved to this scale. 



B. Error Analysis 

All fits in this work are performed by minimizing the 
with respect to the S'-wave event distributions reported 
by the BES PWA [H El [13 for ^tt+tt", wtt+tt- and 
(j)K'^K~. The only parameter left fixed in all of the fits 
is the momentum cutoff (/max, which is constrained by 
the requirement that the /o(980) resonance generated by 
the FSI should be well centered on the corresponding 
peak in the experimental (j)TT~^Tr~ spectrum. A rather 
broad minimum is found around a cutoff of ~ 0.9 GeV, 
and in order to avoid unnecessary fine-tuning, that value 
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TABLE I: Parameter values obtained from a simultaneous fit to the BES results 

Elllllll for J/^ 

and J/ijj ^ Lu-K^-K^ . All values of the are quoted at the scale /i = nip. The meson masses and other input parameters are 
given in Sect. lITr^ and the uncertainties of the parameters correspond to the error analysis described in Sect. lITrSl Fixed 
parameters not determined by the fits are given in parentheses. The fits of type I optimize agreement for J/ijj ^ cj>K^ K~ , 
whereas type II fits do so for J/i/) — > cj)-K'^-K~ . Fits III and IV are similar to "Fit I" in the weighting of the input data. Note 
that "Fit I" represents the main result of this study, as elaborated in Sect. IIIll The dimeson spectra which correspond to the 
main fits (I and II) are shown in Fig.|4] and the corresponding sets of scalar FFs in Figs. and |^ 



Fit LI [10-^] LI [10-^] LI [10-==] LI [10-3] Q [keV"*] g^,, [GeV] {r^Ys [fm^] 



I 


0.84 


+0.06 
-0.05 


0.45 


03+° !^ 

'^•'-'■J -0.13 


S3+°-i* 

^■^^ -0.17 


42.1 


+ 5.0 
-5.0 


0.132 


+0.018 
-0.015 


(0.90 ± 0.025) 


0.657 


la 


0.817 


+0.048 
-0.029 


0.384 


-0 025+°"° 


0.328+0:1?^ 


45.27 


+ 1.82 
-3.07 


0.133 


+0.017 
-0.016 


(0.875) 


0.652 


lb 


0.840 


+0.043 
-0.023 


450 +^-°°^ 


0.028 t°:J?I 


325 +0 129 

U.O^O _0.168 


42.08 


+ 1.24 
-2.49 


0.132 


+0.016 
-0.015 


(0.9) 


0.657 


Ic 


0.858 


+0.037 
-0.022 


520 +° '"''' 

U.OZU _o.oi5 


0.081 ^-.^s 


n ■?9'? +0.134 
U.OZO _o.i66 


39.15 


+ 1.09 
-2.03 


0.131 


+0.016 
-0.015 


(0.925) 


0.663 



II 


n qc,+0 08 


25 +°-2-i 

'-'•^^ -0.29 


06+° " 

'J-'J" -0.08 


0.18 


+0.21 
-0.25 


44.7 


+ 1.2 

-4.4 


104 +°°25 


(0.90 ± 0.025) 


0.659 


III 


94+0 05 

^■^^ -0.04 


40 +° i° 

-0.12 


(0.20) 


0.12 


+0.11 
-0.11 


35.1 


+2.5 
-2.6 


1 90 +0.015 
U.IZO _o.oi3 


(0.90 ± 0.025) 


0.670 


IV 


(0.70) 


55 +° i2 

'J-^'J -0.10 


-0 is+° " 

"•l>^ -0.07 


0.55 


+0.13 
-0.14 


50.9 


+3.8 
-6.4 


162 +°°2* 


(0.90 ± 0.025) 


0.634 



is taken as the central value to be used in all fits. To 
estimate the uncertainty associated with the choice of 
gmax, a range of 50 MeV, between 0.875 GeV 

and (/max = 0.925 GeV, has been considered. Larger 
or smaller values produce an /o(980) peak which is no 
longer well centered on the experimental one in ^tt+tt". 
In order to illustrate the dependence on the choice of 
cutoff (7max, fits for a fixed cutoff, each corresponding to 
the values quoted above, have been given in Table as 
"Fit la - Ic" . 

The statistical errors of the BES data are 
small and have little bearing on the actual sources of un- 
certainties of the fitted parameters. The main sources of 
error accounted for in the present analysis are as follows: 
Firstly, the tr amplitude used in the PWA for wtt+tt- 
published by BES 13] did not take into account chan- 
nel coupling effects, which leads to a considerable un- 
certainty in the 5'-wave contribution close to the KK 
threshold. In view of this, the wtt+tt" data is only fitted 
up to bin 30 and the uncertainty on the parameter val- 
ues introduced by this arbitrary choice is estimated via 
fits up to bins 29 and 31. It is found that the effect of 
this is rather small for most parameters, and all of the 
fits contain a sharp characteristic cutoff in the 5-wave 
TT+TT- distribution close to the KK threshold. Secondly, 
the precise magnitude and shape of the S'-wave contri- 
bution to the conspicuous low-energy enhancement, or a 
peak in the [J7r+7r~ data has to be considered as some- 
what uncertain, as the BES PWA models the a contri- 
bution in terms of different functional forms related to 
a relativistic Breit-Wigner description. Such approaches 
are known to suffer from consistency problems with re- 



spect to the data on tttt phase shifts, and indeed only 
one of the BES fits takes those data simultaneously into 
account. The S'-wave contribution, as determined by the 
BES PWA constitutes about 60% of the total signal in 
the a peak of the wtt+tt" spectrum. The fact that the 
properties of the a peak are strongly influenced by in- 
terference with rescattering effects in the wtt system has 
also been demonstrated in Ref. 0] , where a sequential de- 
cay mechanism with intermediate vector and axial- vector 
mesons, e.g. J/ip — > 6i(1235)7r — > o^tttt was investigated. 

The situation in the BES PWA is similar, as most of 
the remaining events at the a peak are due to the inter- 
ference of the a amplitude with the prominent 61 (1235) 
and p(1450) resonances in the ujt: system. At higher en- 
ergies, the left wing of the strong 13- wave /2(1270) res- 
onance also contributes significantly. The situation for 
the (/)7r+7r~ spectrum is much more straightforward since 
the S'-wave contribution completely dominates the data 
up to the KK threshold, and also because of the much 
weaker interaction in the 07r system. The advantage of 
using the BES PWA is that the direct S-wave contri- 
bution to the fj peak has been disentangled from the 
61(1235) and p(1450) contributions, such that a direct 
fit to the S-wave spectrum is possible without consid- 
eration of strong rescattering effects in the o^tt system, 
which would introduce further uncertainties due to the 
modeling of such effects. The downside is that such un- 
certainties are included instead in the determination of 
the S-wave amplitude in the PWA. 

In view of the abovementioned issues, a conservative 
10% uncertainty in the size of the S-wave contribution 
has been adopted for the wtt+tt" spectrum to estimate 
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FIG. 4: Fit results and BES data for J/tp 4>-k+-k~ , J/tp <t)K+ K~ , J/tp ^ u-k+tv' and J/^j -> ujK+K~ . The datapoints 
are uncorrected for background and detection efficiency, and tiie errorbars represent the statistical error. The direct S-wave 
contribution from the BES PWA and the estimated background are denoted by the dashed and dotted histographs, respectively. 
The solid and dashed curves describe "Fit I" and "Fit 11" to the S-wave contributions. As described in the text, "Fit III" is 
performed for a fixed value of Lg and "Fit IV" for a fixed value of L\. The fitted parameter values are given for each case in 
TableQ] As discussed in Sect. lIIl'ni the J /ip luK'^K~ spectrum is not included in the fits. 



the errors on the fitted parameter values. Moreover, as 
described in Sect. IIII^ the wtt+tt" data has been given 
a rather low weight factor, and thus the shapes of the 
fitted spectra are essentially dictated by the ^tt'^tt^ and 



4>K^K^ channels. Nevertheless, inclusion of the wtt+tt^ 
data in the fits is essential in order to provide informa- 
tion on the relative intensities of the (f)Tr'^n~ and wtt+tt" 
channels. Without that information, the model parame- 
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ter cannot be pro perly constrained since the present 
BES data [H [ll 111 are unnormalized, which was also 
the case in Ref. Q. 



C. Fit Results for (j>Ti+n- and (j>K+ K' 

The experimental dipion spectrum for J /ip (p-K^i:^ 
is dominated by the conspicuous S'-wave /o(980) reso- 
nance, such that an accurate description of that state is 
necessary in order for a reasonable fit to be attained. In 
the present description, the /o(980) is dynamically gen- 
erated and its appearance is a consequence of the channel 
coupling to KK. However, since the FSI model depends 
on only one parameter, namely the momentum cutoff 
9max, the possibilities of optimizing the width of that 
state are severely restricted. It is thus fortunate that 
the present considerations give a reasonable description 
of the /o(980), with the possible objection of that state 
being slightly too narrow. This generates an ambiguity 
which is reflected in the fact that two different fits, shown 
as "Fit I" and "Fit 11" in Fig.^ may be achieved depend- 
ing on whether the (/)7r+7r^ or the (pK^K^ spectrum is 
optimized. 

As in the present case for "Fit 11" , agreement was op- 
timized for (f)TT'^Ti~' in Ref. 0], which resulted in a no- 
ticeable underprediction of the (j)K~^K~ spectrum. As 
explained in Sect. IIII Al "Fit I" was defined by forcing 
agreement with the (j)K^K~ spectrum, which may be 
accomplished at the small price of a moderate overpre- 
diction in the ^tt+tt^ spectrum, which amounts to a fac- 
tor ~ 2 at the bin centered on the /o(980) peak. A com- 
parison of the fitted parameter values given in Table ^ 
suggests that the ones corresponding to "Fit I" are more 
realistic. Further comparison of the two fits has also re- 
vealed that "Fit I" is more stable against variations in 
the input. It is thus concluded that "Fit I" is preferable 
and that the modest overprediction at the /o(980) reso- 
nance in the ^tt+tt" spectrum is due to the limitations 
of the FSI model employed in this paper. It should be 
emphasized that this discrepancy is small enough to go 
unnoticed in a study of meson-meson phase shifts alone. 
Also, the difference between "Fit I" and "Fit 11" is neg- 
ligible in the wtt+tt" spectrum. 

Finally, the observed excess of events in the left wing 
of the /o(980) resonance in the (/)7r+7r~ spectrum should 
be mentioned, which indicates a significant presence of a 
nonstrange, or a contribution. Indeed, the fitted value of 
^ 0.13 given in Table ^ reproduces this feature well. 

D. Fit Results for ijjtt~^tv~ and luK'^ K~ 

The most conspicuous features in tJvr+Tr^ are the lack 
of a significant peak at the position of the /o(980) res- 
onance, and a very prominent a peak around a dipion 
energy of ~ 0.5 GeV. It is evident that "Fit I" in the 
lower left-hand plot in Fig. ^ faithfully reproduces both 



of these features. The a peak is generated by the domi- 
nance of the pion nonstrange FF, which in contrast to the 
strange one approaches a constant value at low energies. 
The main difference between the S'-wave contribution of 
the BES PWA and the present fit is close to the KK 
threshold, where the fitted spectrum is sharply cut off by 
channel-coupling effects, producing a plateaulike appear- 
ance close to the /o(980) resonance. This characteristic 
shape is produced here by the interplay of the strange 
and nonstrange scalar EEs close to the KK threshold. It 
should also be noted that most of the remaining tttt events 
below the KK threshold are, as elaborated in Sect. IIII Bl 
due to interference effects with the 6i(1235) and p(1450) 
resonances in the ujtt system. 

Once the parameters and are fixed by a fit to 

the 07r+7r~, wtt+tt" and (j)K^K~ spectra, the scalar EEs 
can be used to predict the S'-wave contribution in the 
LuK^K~ channel. It turns out that such a prediction 
is quite robust, since attempts at tweaking the fit by 
forcing changes in the u)K^K~ sector do not have any 
appreciable effect. Nevertheless, it is immediately appar- 
ent from Fig. ^ that the S'-wave contribution reported 
by the BES collaboration does not compare favorably 
with the present results. The S'-wave fitted by BES con- 
tributes ~ 50 events close to the KK threshold, whereas 
the present prediction can account for ~ 20 events. This 
number can be motivated from the present considera- 
tions as follows: From Figs. |S1 and El the moduli of the 
nonstrange pion and kaon EEs are seen to be approxi- 
mately equal in magnitude, such that a roughly similar 
S-wave contribution would be expected at ~ 1.1 GeV 
in the wtt+tt" and ujK~^K~ spectra. This would lead 
one to expect ~ 50 events in the ujK~^K~ spectrum as 
well, however the detection efficiency reported by BES for 
LoK^K~ is only about half of that reported for o^tt+vt^, 
yielding the calculated ~ 20 events. The situation in 
the (j)K^K^ spectrum is different, as the modulus of the 
strange kaon FF at ~ 1.1 GeV is much larger than that 
of the strange pion EE, as evidenced by Figs. |S1 and El 

With respect to this discrepancy, it should be noted 
that when the BES PWA was done, no reliable infor- 
mation was available on the magnitude of a possible a 
coupling to KK, and therefore it was fitted freely in that 
analysis. The BES data exhibits interference between 
a KK and a signal fitted as iiri(1945) loK. It 
turns out that the latter can be adjusted slightly in mass 
and width in order to reduce the a — *■ KK signal to a 
level consistent, within statistical and systematical er- 
rors, with the findings of the analysis presented in this 
paper, at the price of a slightly worse fit to the KK 
event distribution^. The overall changes induced in the 
BES PWA by such modifications are very small. 



^ The authors thank David Bugg for his kind cooperation in inves- 
tigating the discrepancy between the predicted u)K^ spec- 
trum and the BES analysis. 
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FIG. 5: The nonstrange and strange scalar FFs F" and F^ of the pion. Moduh, real and imaginary parts are shown for the 
parameters corresponding to "Fit I" and "Fit H" in Table Q The curves labeled "Old" show, for comparison, the scalar FFs 
for the parameters obtained in the analysis of Ref. (3- 



E. Analysis of the L[ 

One of the main objectives of this analysis is the de- 
termination of the LECs L4 and Lg, for which the data 
+ PWA provided by the BES collaboration is certainly 
precise enough. As discussed in Refs. Tj','35'|, where a set 
of standard estimates for the LECs of CHPT are given. 



L4 and Lg are expected, in terms of large Nc arguments, 
to vanish at some unknown scale in the resonance region 
of QCD, conventionally taken to be nip or m^. As dis- 
cussed above, "Fit I", as given in Table is considered 
the main result in this respect. The values so obtained, 
(0.84 ± 0.06) X 10-3 for LI and (0.03 ± 0.16) x 10-^ for 
Lg should be compared with the NNLO CHPT analy- 
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ses of Refs. "33*, "3^ , where the L\ were extracted by 
fits to the available experimental data on Ki4^ decays. In 
that study, constraints on L\ and ig were also derived 
by requiring a properly convergent behavior of the chi- 
ral expansion up to NNLO for several quantities, such as 
the pseudoscalar meson masses and decay constants. It 
was found that such a constraint can only be satisfied in 
a rather small region centered around L\ — 0.2 x lO^'^ 
and Lg = 0.5 x 10"'^. Also of interest in this context 
is the analysis based on QCD sum rules in Ref. [3^ . 
where Lg was constrained to be positive and in the region 
0.2 X 10-3 < < 0.6 X 10-3. Also, a significantly P2S- 
itive value of L\ = 0.4 x 10-^ was obtained in Ref. |37|. 
although with unknown error. Some updated results for 
the L\ are given in Ref. {i^l , and a recent analysis of -kit 
and TiK scattering is presented in Ref. :39j, where the 
preferred values of L\ and ig were found to be compat- 
ible with zero. Further recent determinations of the L\ 
include the analyses of Refs. ^29, 30] and 4^ in terms 
of the Inverse Amplitude Method (lAM) of Ref^ lij . In 
Ref. a simultaneous fit to the tttt and KK partial 
wave amplitudes for / = 0, 1, 2 was performed using the 
complete NLO CHPT amplitudes, which yielded a value 
oiLl equal to (0.2 ±0.1) x lO^^. The analysis of Ref. [13 
considered all two-meson scattering amplitudes includ- 
ing the 7777 channel within the chiral lAM framework, 
which enabled the simultaneous extraction of the LECs 
L\ through Lg, the reported values of L\ and Lg being 
(-0.36 ±0.17) X 10-3 and (0.07 ±0.08) x 10-^, respec- 
tively. Finally, Lattice QCD also offers the possibility of 
determining the LECs of CHPT via application of Par- 
tially Quenched (PQ) Lattice QCD and PQCHPT. For 
a theoretical background, see Refs. [H ElEl El Ell, 
and some recent Lattice QCD results for the L[ can be 
found in Refs. [Uliilil. 



F. Results for L\ 

The value of L\ preferred by the present study is 
(0.84 ± 0.06) X 10^3^ which is largely consistent with the 
value (0.71±0.11) x 10^3 reported in Ref. 7] where a sim- 
ilar analysis was performed, although with much poorer 
data and without the guidance of a comprehensive PWA. 
These results are about twice as large as the largest val- 
ues found by the other analyses presented above. Also, 
the present result for i2 is difficult to reconcile with the 
lAM analysis of Ref. [33 , where a rather large but neg- 
ative value of L4 ~ —0.4 x 10-3 was obtained. However, 
the results of that lAM analysis are not very sensitive to 
the value of L\, since it has been pointed out in Ref. [sJl 
that values of L4 ~ or L\ ~ 0.2 x 10-3 can be accom- 
modated with small overall changes. 

The present value of L\ is constrained by a simultane- 
ous fit to the (/fTT+TT- and uj'jT^'n~ spectra. Even a modest 
deviation from the central value of L\ tends to upset the 
interference between the strange and nonstrange FFs at 
the /o(980) peak, producing for large deviations a pro- 



nounced cusp in the wtt+tt- spectrum which contradicts 
the flat, plateaulike behavior seen in the experimental 
data. Nevertheless, the wtt+tt^ spectrum is not patholog- 
ically sensitive to the exact properties of the /o(980) res- 
onance, since the bins centered on that energy in tJTr+Tr^ 
have been excluded from the fit; Rather, a smaller value 
of L\ as used for "Fit IV" , starts to upset the behavior 
of the low-energy parts of both the (j)'!T~^'!T~ and wtt+tt" 
spectra. It is thus encouraging to note that once is al- 
lowed to adjust itself optimally such that the low-energy 
region is well described, then the o^tt+tt^ spectrum also 
"finds" the right shape at the K^K^ threshold. The 
conclusion is therefore that a smaller value of L\ can- 
not easily be accommodated within the present analysis. 
Since the OZI rule is also based on large Nc arguments, 
a large value of L\ may be interpreted as a clear signal 
of OZI violation. However, it should be emphasized here 
that the inclusion of higher-order terms in the vertex of 
Eq. ^ may affect the value of L\ significantly, although 
as discussed in Sect. IIII Jl there is little evidence as yet 
for such higher-order contributions. 



G. Results for Lg 

For LI, the value (0.03 ±0.16) x 10-3 determined from 
"Fit I" is in much better agreement with the results of 
the other analyses above. In remarkable contrast to the 
situation for L\, the present value of Lg is completely 
consistent with that determined by the lAM method in 
Ref. [331 . Although the present result is consistent with 
zero within errors, a fixed value of around 0.2 x 10-3, as 
suggested by the bounds on Lg derived in Ref. [sj , does 
little to upset the agreement with experiment, as indi- 
cated by "Fit HI", the main effect being a slight over- 
prediction close to the /o(980) resonance in the c^tt+tt- 
spectrum. 

As already noted in Sect. IIIl'Hl the /o(980) as gener- 
ated by the FSI is slightly too pronounced in the present 
framework, which may be interpreted as pulling the fit 
toward lower values of Lg. Anyway, as the main sensi- 
tivity to Lg lies in the region of the /o(980) resonance, it 
is unfortunate that the presently available version of the 
BES PWA does not take into account channel coupling 
effects for that state. If the bins of wtt+tt- data around 
the /o(980) could be included in the fit, the uncertainty 
on Lg would most likely be much reduced. The value of 
Lq found by the old analysis in Ref. 018-0.22x10-3^ 
without any error specified, which is inconsistent with 
the lower bound given in Ref. js^l- However, the ba- 
sic reason such a value was obtained is that the wtt+tt- 
spectrum was assumed, due to the lack of a PWA, to 
be completely dominated by the S'-wave contribution at 
low energies. However, this may only constitute about 
~ 60% of the total event count, as indicated by the re- 
cent BES PWA. For this reason, the wtt+tt- spectrum as 
given in Ref. 7] does not have the proper shape close to 
the KK threshold. 
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FIG. 6: Tlie nonstrange and strange scalar FFs and r| of the Icaon. Moduli, real and imaginary parts are shown for the 
parameters corresponding to "Fit I" and "Fit 11" in Table |I] The curves labeled "Old" show, for comparison, the FFs for the 
parameters obtained in the analysis of Ref. 



H. Results for L\ and L\ 

The values of L\ and L\ are much less controversial, 
as most studies are in reasonable agreement concerning 
their numerical values. In the previous work of Ref. Q 
their values were kept fixed and were taken to coincide 
with those determined by the Ki^ analysis of Ref. [s^. 



These values were quoted as L\ — (0.65 ± 0.12) x 10"^ 
and Lg = (0.48 ±0.18) x 10~^, respectively. However, in 
the present analysis it turns out to be possible to let these 
parameters adjust themselves to their optimal values. 
From "Fit I" presented in Tabled the obtained values are 
LI = (0.45 ± 0.09) X 10-3 and = (0.33 ± 0.17) x lO'^. 
Although they come out slightly smaller than the values 
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of Ref. [33 , it is encouraging to note that an optimal fit to 
data does not require extreme values of Lg and Lg. The 
more recent fit in Ref. [s^l has reported an updated set 
of LECs which compare slightly more unfavorably with 
the present values, namely = (0.91 ±0.15) x 10"'^ and 
LI = (0.62 ±0.20) X 10~3. However, it should be kept in 
mind that the nature of the fits in Ref. "s^ did not allow 
for a detailed consideration of the effects of the LECs 
and Lg. It should also be noted that the values for L5 and 
Lg obtained in Refs. js^ls^ are significantly smaller, and 
therefore in much better agreement with the present re- 
sults, than the old standard values = (1.4±0.5) x 10~^ 
and LI = (0.9 ±0.3) x lO'^ in Ref. [H. Finally, it should 
be noted that the results of the lAM analysis of Ref. j^O] 
for Lg and Lg are essentially identical to the standard 
ones given in Ref. |3^. 



I. Fitted Values for and 

The overall normalization constant is determined 
by the present analysis to be around 40 keV^^^ as in- 
dicated in Table HJ This value is difficult to compare 
with Ref. 7], as that work did not consider the detec- 
tion efhciencies of the DM2 and MARK-III experiments. 
Furthermore, Ref. also employed a different definition 
of the overall normalization factor. In the present study, 
the dimension of is a consequence of the absorption 
of the proportionality factor in Eq. 1)15(1 into the product 
V2 gBoC^. This factor could be estimated once a con- 
version of the event distributions reported by the BES 
collaboration to differential branching fractions is avail- 
able. However, such an estimate is complicated by the 
nontriviality of the detection efficiencies shown in Fig. 13 
and by a considerable uncertainty in the hitherto pub- 
lished values for the relevant branching fractions [Sj . In 
principle, could then be used, as in Ref. to extract 
a value for the dimensionless product gBo. A determi- 
nation of the coupling constant g itself is more difhcult 
since the precise value of Bq , which is related to the quark 
condensate, is only approximately known. In view of the 
above issues, further analysis of is beyond the scope 
of this paper. 

Of considerably more interest is the fitted value of the 
parameter A^, which directly measures the OZI violation 
in the J/4> ^ (/itt+tt" and related decays. From "Fit I" 
in TableQI the present value is A^ — 0.13 ±0.02, which is 
lower than the value 0.17 ± 0.06 obtained in Ref. ^, al- 
though both values are consistent given the rather large 
errors of the latter result. The reason for this apparent 
decrease in A^ can be immediately understood in terms 
of Eq. ^ , if one notes that all of the events in the low- 
energy part of the unr'^7r~ spectrum were taken in Ref. 
to be iS-wave events. Put another way, the 5'-wave con- 
tribution for o^TT+TT^ was taken to be about 50% larger in 
Ref. 't'I than it was found to be in the BES PWA. A larger 
value of was thus required, in turn yielding a larger 



value of A^ via Eq. (jSj). Nevertheless, the present value of 
A^ indicates a significant amount of OZI violation, which 
manifests itself as a contribution to (/)7r+7r~ from the non- 
strange scalar FF, and a contribution to wtt+tt" from the 
strange one. In particular, a realistic description of the 
WTT+TT" spectrum close to the KK threshold requires an 
interplay between the nonstrange and strange scalar FFs. 

In this context, it is useful to recall Ref. Q, where 
the experimental results of Refs. l9l lTollTl| were normal- 
ized using the available information l34| on the width 
and branching fractions of the J/i(j. Since Ref. con- 
sidered only the LO CHPT contributions to the oper- 
ators fin and ss, the values of the do not appear 
as additional parameters in that model. It should be 
pointed out though, that the FSI between pairs of pi- 
ous and kaons was taken into account in a similar way 
as in Ref. and in the present work. Furthermore, the 
strength of the sequential decay mechanism via interme- 
diate vector and axial-vector resonances could be fixed 
from experimental information js^l up to a relative sign 
of the respective coupling constants. Up to this rela- 
tive sign, the only free parameters in the calculation of 
Ref. 13 were gBo and A^. The sign ambiguity could not 
be resolved in Ref. , since equally good fits to the data 
were obtained for both choices. If both the J/ip VP 
and J/ijj AP coupling constants have the same sign, 
Ref. la finds gB^ = 0.032 ± 0.001 and A^ = 0.12 ± 0.03. 
This value of A^ is consistent with the one obtained from 
"Fit I" in Table n 



J. Results for the Scalar Form Factors 

Once all parameters are determined by the analysis of 
the BES data 0, 0, 0|, the shapes and magnitudes 
of the scalar FFs of the pion and the kaon are fixed. 
In Figs. El and [HI they are shown for the fitted param- 
eters given in Table P and compared to the form fac- 
tors obtained in Ref. 7]. It is evident that the values of 
the form factors at the origin have changed significantly 
with respect to those of Ref. i7|, even though the gross 
features remain similar. The most significant change is 
that the form factors now have a realistic behavior close 
to the KK threshold, which is reflected in the values of 
the form factors at the origin. Thus r"(0) has increased 
by ^ 30%, while the value of Ff (0) remains very small 
and negative. The situation for the kaon form factors is 
more drastic, since they were not very well determined in 
Ref. 7]. The present results indicate that r5(0) is small 
and negative, while r2(0) is positive and has increased 
significantly from the value determined by previous work. 

Since the description of the experimental data shown in 
Fig.2|is satisfactory, a natural question to ask is whether 
the scalar FFs can be considered realistic as well. The 
issue is complicated by the essentially unknown momen- 
tum dependence of the vertex given in Eq. , for which 
the simplest possible form was assumed in Ref. 7j. The 
question is then whether the scalar FFs are simulating 
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some neglected momentum-dependence in the J/tp ^ (j) 
and J/tp ^ ui vertices. Although the detailed study of 
this is beyond the scope of the present study, the scalar 
radius (r^)^ of the pion and the curvature of the non- 
strange scalar FF may nevertheless be used to check for 
consistency with other theoretical descriptions. These 
quantities are defined |l8l l50j | in terms of the expansion 

^ = l + -(r2)^. + <.2 + 0{s'), (42) 

of the nonstrange scalar FF of the pion around s — 0. 
The values for (r^)J given by "Fit I" in Tableware quite 
stable around 0.657 ± 0.006 fm^, whereas the values for 
Cs are - 12.50 GeV"'' for "Fit I" and ~ 12.53 GeV* 
for "Fit 11". These values can be compared with other 
theoretical predictions, starting from the prediction of 
Ref. J.7J, where an estimate of (r^)J = 0.55 ± 0.15 fm^ 
was obtained from the experimental value of Fk/ F-^, ne- 
glecting OZI-violating contributions. This estimate has 
been refined by the dispersive work of Ref. |5lj| and fur- 
ther improved upon in Ref. 52] using the tttt Roy equa- 
tions to yield (r^)^ = 0.61 ± 0.04 fm^. In Ref. [p, the 
curvature was also determined and found to be around 
^ 10.6 GeV"'', which is about 10% smaller than the 
present value. In the NNLO CHPT analysis of Ref. pl |. 
the value of could not be pinned down very pre- 
cisely, although the preferred values are in the range 
10... 13 GcV-4. 

The effects of inelasticities due to the tttt Att and 
TTTT rjrj channels were studied in detail in Ref. [3^ . 
which produced a ran ge o f values for (r^)^ consistent 
with the result of Ref. |52| . Higher-order CHPT correc- 
tions to the scalar FFs have been considered in Ref. |50l | 
and in the NNLO fits of Refs. js^jH^I, which also yielded 
a preferred value of (r^)^ around ~ 0.61 fm^. Finally, the 
value of {r'^Ys ^-l^o been extracted from the S'-wave 
TTTT phase shift for / = using the sum rule of Ref. j 
which yielded a value of 0.75 ± 0.07 fm^. However, as 
stated in Ref. 54], this work is not in agreement with 
earlier determinations of {r'^)1- 

The conclusion of this comparison is that the present 
values for the pion scalar radius and the curvature 

of the pion scalar FF are largely consistent with other 
recent theoretical descriptions. The present (r^)J value 
of ^ 0.65 fm^ is at the upper limit of the generally ac- 
cepted value of 0.61 ± 0.04 fm^. It is possible that the 
slightly higher value found in the present study is hint- 
ing at a neglected momentum dependence in the vertex 
of Eq. ^ or at some inaccuracy in the present descrip- 
tion of the TTTT phase shifts. Nevertheless, the reasonable 
value of is reassuring in view of the somewhat un- 

orthodox value of L\ obtained from the present analysis. 

IV. DISCUSSION AND OUTLOOK 

The main results of this analysis of the experimen- 
tal EES results [H El H are the set of values ob- 



tained for the L[ and the corresponding sets of scalar FFs 
of the pions and kaons. Due to the accuracy of the 
BES data, the actual errors of the are dominated 
by the uncertainty in the theoretical description. The 
present value of Lg is small and compares well with the 
known theoretical bounds. The present value is com- 
patible with zero, although on the positive side. On the 
other hand, L\ was found to be large and positive, around 
(0.84±0.06) X 10~^. It has also been shown that a value of 
L\ smaller than about 0.7x 10~^ is not compatible with 
the present analysis. Although this value of L\ is in con- 
flict with most of the presently available fits of the L[, an 
advantage of the present analysis is that the expressions 
for the scalar FFs depend on only four LECs, namely 
L4,Lg,ivg and Lg, and that no assumptions about their 
values had to be made, since all of them could be con- 
strained in a unique way. Also, it should be kept in mind 
that the present description of the scalar FFs, although 
certainly not complete, contains some of the features that 
would be present in an all-orders calculation in CHPT. In 
contrast, determinations of the L\ from truncated NLO 
or NNLO CHPT expansions suffer from possibly large 
errors due to the neglect of higher-order contributions, 
and the poorly known values of the NNLO LECs. 

Naturally, a further question to ask is whether the fit- 
ted value of L\ could be modified by changing some of 
the features of the present description of the J/^ decays; 
Possible avenues of further investigation include the con- 
sideration of the "qri channel in the FSI which was ne- 
glected in this work, although it cannot be argued that 
this is necessarily a large effect as the description of the 
FSI is independent of the L\ at this level of theoretical 
sophistication. It should also be noted that in the lAM 
framework of Ref. 41], the tttt phase shift is given in 
terms of the and thus more fiexibility of adjustment 
is possible. On the other hand, such an approach may ac- 
tually turn out to be more constraining than the present 
one, since the pole position of the /o(980) resonance is 
then directly affected by the L\. 

Another essentially open question is whether a more 
relevant decay vertex than the one given in Eq. should 
be considered or not. In order to give a partial answer 
to this question, the pion scalar radius has been esti- 
mated for the present set of scalar FFs. If the FFs were 
compensating heavily for some neglected momentum de- 
pendence in the J — > 0, w vertex, this might reveal 
itself as an unrealistic value for The present anal- 

ysis gives a value close to ~ 0.65 fm^ which, as argued in 
Sect, nil Jl is in reasonable agreement with other recent 
theoretical predictions. It therefore appears that there is 
no compelling evidence as yet for a significant momentum 
dependence in the vertex described here by Eq. (^) . 

An issue often raised is the different appearance of the 
dipion spectra in the ldtttt and ^tttt channels, which are 
shown in Fig. 0| Within the present framework, this dif- 
ference can be understood in very simple terms. Accord- 
ing to the expectations from the OZI rule, the wtttt spec- 
trum should be dominated by the nonstrange scalar FF 
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of the pion, whereas the ipTTir process should essentially 
be driven by the strange scalar FF alone. However, the 
latter one vanishes to LO in CHPT, which leads to a 
strong suppression of the low-energy tail of the dipion 
spectrum in (j)TnT. On the other hand, the nonstrange FF 
is nonzero to LO in CHPT, and thus the decay amplitude 
itself does not vanish for low dipion energies in the ujittt 
decay. The appearance of the experimental lutttt spec- 
trum near the tttt threshold is then given by the squared 
modulus of the nonstrange scalar FF and phase space. 

In this context the OZI violation parameter is in- 
structive, since it determines whether any contribution 
from the nonstrange pion scalar FF is present in the (/)7r7r 
decay process. It turned out that could be quite accu- 
rately determined from the present fits, which yielded a 
value of A^ = 0.13 ± 0.02, which indicates that the non- 
strange pion scalar FF contributes significantly to the 
(jyTTTT decay, and indeed this contribution appears as a 
distinct enhancement in the low-energy tail of the di- 
pion distribution shown in Fig. 21 Nevertheless, this 
parameter has decreased from the old value of around 
A^ ~ 0.17 given in Ref. 0, due to the fact that the S- 
wave (T contribution to the lotttt decay has turned out to 
be smaller than previously estimated. As discussed in 
Sect. llllll the present value of A^ is consistent with that 
obtained in Ref. where the strong interaction in the 
ujir and (/)7r systems was explicitly taken into account. 
This was not done in the present work since a deter- 
mination of the direct S'-wave contribution via PWA is 
available _12, 13, 14j. 

An extra bonus of the present work is that the dikaon 



distribution in the ujK^K~ decay could be predicted 
close to the KK threshold, which is instructive since the 
present experimental data for that channel is difhcult to 
interpret in terms of a PWA Ti^. In the same way, the 
predicted shape of the S'-wave distribution in the uj'it'^-k~ 
decay is slightly different from the one favored by the 
BES PWA. Finally, as the scalar FFs of the pion and the 
kaon have been strongly constrained by the J decays 
analyzed in this paper, the impact of these newly deter- 
mined scalar FFs can be studied for other processes like 
B ^ (TTT [s^ or B decays into three pseudoscalars, such 
as B ^ Kg-mr [H^. Work along these lines is under way. 
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